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O ■ Abstract 

■ We investigate the azimuthal distribution of quark-antiquark jets in DIS difFractive dissociation 
with large transverse momentum. In this kinematical region the matrix element is expressed 

^ • in terms of the gluon structure function. For the transverse part of the cross section we find 

I I a cos(20)-distribution with the maximum at (/> = ±7r/2, i.e. the jets prefer a direction perpen- 

dicular to the electron plane. This is in contrast to boson gluon fusion where the qq jet cross 
(-H I section for transversely polarized bosons peaks at = and (j) ^ n. We discuss the origin of 

this striking difference and present numerical results relevant for the diffractive dissociation at 
HERA. 

X' 
J3 



1.) The observation and analysis of diffractive events in deep inelastic electron proton scattering 
at HERA ^, ^, ^ has attracted much attention recently. Particular interest has been given to 
the question whether this process is dominated by the same (soft) Pomeron which is observed in 
difFractive dissociation of hadron-hadron scattering, or by the hard Pomeron which is associated 
with the observed rise of F2 at small x. As far as the (inclusive) diffractive structure function is 
concerned, one expects to see a combination of both the soft and the hard Pomeron. As a method 
of separating the hard component it has recently been suggested to look for final states of the 
difFractive system which consist of jets with large transverse momenta. The simplest configuration 
is a two-jet final state of a quark-antiquark pair. The cross section for this process, together with a 
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few numerical predictions, has been presented in |5|. In this letter we continue this analysis by cal- 
culating, in the 7* — Pomeron center of mass system, the azimuthal distribution of the jets relative 
to the plane formed by the beam axis and the scattered electron. As we will show, this distribution 
differs from that of the photon-gluon fusion process and therefore presents a characteristic signal 
of the two gluon model for the hard Pomeron. 



2.) The process that we are going to investigate is shown in Fig. la, and the variables are il 



lustrated in Fig. lb; for simplicity we restrict ourselves to the forward direction t = 0. Throughout 
this paper we work in the Pomeron-photon center of mass system. We begin with the subprocess 
7* + proton qq + proton, and we take the incoming photon to be linearly polarized in the trans- 
verse direction. The polarization vectors are = (0, 1,0,0) and ey = (0,0, 1,0), and the electron 
momenta lie in the x-direction. We work in the leading-log (1/x) approximation and use the k- 
factorization theorem Q to express the amplitude through the unintegrated gluon distribution of 
the proton. This is based on the assumption that the total energy s = {q + p)^ is much larger 
than the photon virtuality = —q^ and the invariant mass of the outgoing quark anti-quark 
pair = {q + xpp)^ . The quark phase space is parameterized in terms of the quarks transverse 
momentum k±_ = \k_\_\ {cx cos (j) + Cy sin (j)) , and its longitudinal momentum fraction a according to 
the Sudakov decomposition A;^ = aq'^ + k^/(2apg')p^ + fc^^ with q' = q + xp and = —k±^k'^. 
The invariant mass M"^ is related to the transverse momentum through k'^ = a(l — a)M^. The 
differential cross section for the incoming photon with a polarization lying in the electron plane has 
the form: 

- "^D'T m 26 (I) 

dM^dk^dtd(l)\t=Q dM^dk^dt\t=o 1-2^ ^ dM^dk^dt\t=o 

* 

^^^^^ dM'^'dk^dt ^^^^ given below in eq. (P). The main feature to be noticed is the minus sign 
in front of the angular dependent piece of the cross section. It is opposite to the sign structure of 
the photon-gluon fusion process 0] (which will be discussecd below), and it leads to a character- 
istic maximum of the cross section^ at 4> = tt/2. For the incoming photon being polarized in the 
^/-direction the cross section has the same form as (1), but with a plus sign in front of the second 
term. 



Turning to electroproduction, we have to sum over all polarizations of the virtual photon; the 
interference of longitudinal and transverse polarizations leads to an additional azimuthal depen- 
dence. We need the Sudakov parametrization of the electron momentum: 



V ^ y 



(2) 



Inserting this into the usual lepton tensor: 



Lp,p — 2 



(3) 



similar phenomenon has been observed in the azimuthal distribution of forward jets in deep inelastic electron 
proton scattering ^ . 
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and contracting with the photon polarization vectors we arrive at the following flux factors for the 
linearly polarized transverse photons 



2 ^ 



Summing over all polarizations we arrive at the ep-cross section: 



Op 



dA/Pdk'^dt\t=o 



2(1 -y) 



k2 

IP' 



1 



cos 2(j) 



(4) 
(5) 



daj) J, 
dMhMdt\t-- 



dM^dk^dt \t=o ^ ^dM2(ik2dt it=o 



where the indices T, L, and I refer to the contributions of transverse and longitudinal photons and 
the interference term, resp. Apart from the flux factors of the photons, the transverse part has the 
same structure as seen in eq. (|l]). In particular, the minus sign in front of the angular dependent 
term is a direct consequence of combining (|^) with the flux factor in (Q) . 

Eq.(6) is written for the case in which the polar-angle 6 (cos 9 = 1 — 2a), the angle between 
the quark jet with momentum k and the proton (Pomeron), is restricted to be smaller than '7r/2, 
i.e. a varies between and 1/2. If this jet lies in the other hemisphere (polar angle 9 between 7r/2 
and TT or, equivalently, a between 1/2 and 1), the last term in (6) changes sign. When averaging 
over both hemispheres the cross section above has to be multiplied with a factor of two, and the 
interference term disappears. Then, the only angular dependence results from the term propor- 
tional to cos2i;^) in eq. (^). Note the additional factor in front of this term. We therefore expect 
the azimuthal asymmetry to be more pronounced at large transverse momenta (at the same time, 
however, the cross section decreases with growing transverse momentum). Finally, by integrating 
the cross section from to 27r we recover our previous result |^ . A numerical analysis of the above 
formulae will be performed in part 4.). 

The expressions for the longitudinal and transverse photon-proton cross sections have been ob- 
tained in 101: 



27r 



27r 



1 k^ 
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dM^dk^dt\t=o M^Q^ /i_4k2"3^'^'"'"^''^"' 
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dajjj, 
dWdkMt\t=o 
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lT{Q'',M^y 



(7) 
(8) 



In the same spirit we find for the interference contribution 



27r 



,2 \ 2 



dM^dk^dt\t=o M^KQ"^ 6 



174 7^ -5]e^aemVr2«.'/r(Q',M2,k2)/i(Q^M2,k2) 



(9) 



with 



It 



II 
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12 



+ 



12 + ^(Q2_M2 



M2-Q2 



k2Q2 



and the unintegrated gluon distribution of the proton 

j dX^rdxiP, 1') = xpG{xjp, Q2) 



•^G(a:jP,l2) 
^G(xiP,l') 



IP 



(10) 
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(12) 



In our expression (|T^) we have used, for the two-gluon subamphtude, the (forward) gluon struc- 
ture function of the proton. Strictly speaking, this is not quite accurate, because the longitudinal 
momenta of the two gluon lines differ by xpp. In our leading-log calculation, however, we cannot 
distinguish between this 'slightly nonforward' gluon structure function and the usual DIS gluon 
structure function. We therefore shall use the usual gluon density of the proton, and, as a con- 
sequence, we have to accept an error in the absolute normalization which is characteristic for a 
leading-log calculation (cf. the discussion in |^). 

Without going into more detail we quote the results that were obtained in Q for the functions 
It, Il above 



It 
II 



+ bT- 



d 



k2(M2 + Q2)3 ^ Ql^2 
Q2m2(Q2 _^2) 



xjpG(a;jp,k' 



2Q2 + M2 



k2(M2 + Q2) 
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xipG{x]p,'k' 



M2 ' 

,Q^+M2 
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(13) 
(14) 



In these expressions the first terms represent the double leading log approximation. As discussed 
in the derivative terms with the functions 6t) of and M"^ are those next-to- leading order 
(in logk2(M2 -|- Q'^)/M'^) corrections which, as we believe, are numerically most important. They 
can be found in eqs. (12), (13) of and will not be given explicitly here. We wish, however, 
to stress that these next-to-leading order corrections are not complete, and we believe that a sys- 
tematic study of order-a^ corrections to our 'Born approximation' remains an important future task. 



3.) It will be useful to compare our results with the usual photon gluon fusion cross section. 
To be definite, let us, again, consider the cross section for the production of a quark-antiquark pair 
in the kinematic region where 1^2 much larger than . The momenta are labelled as shown in 
Fig. ^ The gluon momentum has the Sudakov decomposition I = r]p+l± with r/ = xb(1+W^^/Q^) 
(the component along the momentum q' is small and can be neglected), and the polarization vector 
of the gluon is Choosing the incoming photon to have the transverse polarization in 

the x-direction, and summing over the quark helicities, we obtain for the square of the subprocess 
7* + 5 ^ QQ- 

• V • V - 4a(l - a)e^ • V e^. • V (15) 
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where we have disregarded overaU constants, and the vector V is defined as 



As can be seen from (16), the vector V can be interpreted as the qq component of the wave function 
of the transverse photon. It has the form of a dipole. The standard result for the photon-gluon 
fusion process is obtained by integrating over the azimuthal angle (pi of the gluon and taking the 
limit 1^ — > 0. For small P one finds 




M^[D{k)]^ I ^^^^ 



and the cross section becomes 



The angular part now has a positive sign, quite in contrast to the previous case ((|1|) or (|6D). This 
difference can be traced back to the 6ij piece in eq. (p!7[): the second term ~ kjkj alone would have 
lead to the same structure as (|l|), and it is the first term which leads to the sign change in (p!8|). 
As an attempt to find a physical interpretation, one might interpret the contribution ~ kjkj as an 
incoherent product of the two dipoles, whereas the term ~ 6ij represents a correlation between the 
dipole and its conjugate which results from the angular integral. 



For comparison we write down the analogue of eq. dig) for the diffractive case, i.e. the expres- 
sion for the square of the subprocess 7* + 2g — > qq (Fig. |^). It has the same form as (|l5|) , but 
now we have two 1-integrals (both for the matrix element and its complex conjugate). The wave 
functions V in (jl^) have to be replaced by V/)(l) and V£)(l'): 



eo. • Vd(1) ■ Vb(I') - 4a(l - a) • Vb(1) • Vz)(l') (19) 



where 



The limit P — > and the angular integration are done independently for both 1 and 1'. For the 
leading term in this limit we find 

As a result, the tensor VjVj is proportional to kjk,-, and there is no contribution proportional to 6ij. 



To complete our review of the photon-gluon fusion process, we mention that - unlike in the usual 
discussion where only the leading term of the limit 1^ — > is retained - a consistent application of 
the ky-factorization allows to consider also the case F 7^ 0. As one can see immediately from the 
kinematics in Fig. in this case the momentum transfer from the proton system is no longer zero. 
When integrating over 1, one at the same time averages over the transverse momentum of one of 



5 



the two jets, and as a result the cross section can no longer be compared with the difFractive cross 
section of the first part. Nevertheless, we still can ask for the azimuthal dependence of the jet with 
fixed momentum k^. Having performed the 4>i integral, one arrives at a tensor structure similar 
to (0), but the coefficient functions of the two terms might be quite different. Consequently, also 
the sign structure in the cross section in (p^), which results from a combination of these coefficient 
functions, may change, and the azimuthal dependence of the cross section will depend upon the 
1-dependence of the gluon distribution. A definite prediction could be made if one uses the BFKL 
Pomeron as a model for the gluon structure function. In analogy with the results of |^ one expects, 
as a correction to the DLA result ([l8t), a term containing the unintegrated gluon distribution. De- 
tails of such a calculation will be presented elsewhere |I0|. 



4. In this section we want to present some numerical results based on the formula eq. (P). In 
our analysis we closely follow the strategy described in ||5|. For the diffractive cross section the 
integration over t is performed after multiplication of our t = expression with the Donnachie- 
Landshoff formfactor. We use the GRV next to leading order parameterization of the gluon 
density and take account of the subleading corrections, indicated as 6t,6l in eqs. (|T3|), ([l4|). For 
the nondiffractive case where we do not calculate corrections of this kind, we take the GRV leading 
order density. In both cases we take Og running with the scale equal to the momentum scale of the 
gluon density. 

First, to get an overall impression we show in Fig. ^ as a function of the azimuthal angle cp, the to- 
tally integrated electron proton cross section with the kinematical constraints > 10GeV^,xjp < 
10-2, k2 > 2GeV2 and SOGeV^ < W < 220 GeV^. Note that, for giv en values of and we 
have summed over the two configurations with a = ^ + ^ ^/T^^'^^JM^ and a = ^ — ^ ^/l — AH^/M'-^, 
i. e. over the two different hemispheres. In addition we added the degenerate contribution with 
(j) + IT, i. e. the total cross section is recovered by integrating da/d(j) from (p = Q to (j) = tt. 
As a result, the interfence term in (6) drops out, and our curves are symmetric w.r.t. (j) = tt/2. 
For comparison we also show the contribution of the longitudinal polarization (dashed line) and 
the independent part of the transverse polarization (dotted line). The angle dependent contri- 
bution leads to a factor of approximately 2 between the minimum of the cross section at (/> = 0, vr 
and the maximum at = 7r/2. Furthermore, one notices that the longitudinal polarization gives 
on the whole a 10% contribution to the total cross section. As to the differential cross section, 
the relative magnitude of the maximum, as well as the contribution of the longitudinal polariza- 
tion, will depend on the kinematic variables, especially j3 = Q"^ / {Q"^ + M"^) (cf.the discussion in Q). 

Next we want to exhibit the influence of the kinematical prefactors on the relative magnitude 
of the azimuthally asymmetric contribution. In Fig. ^ we show again the totally integrated ep - 
cross section with the above given cuts for Q'^,X]p and W but with three different cuts in k^. The 
curves demonstrate that a restriction of the phase space to larger k^ leads to an increase of the 
asymmetry, due to an additional k^ in the prefactor of the angular dependent term. The absolute 
magnitude of the cross section, however, decreases due to the overall suppression at large transverse 
momenta. 

Instead of k^, it might be more convenient to use the polar angle 9 (Fig. 2) between the jet 
with momentum k and the Pomeron momentum in the photon-Pomeron cm-system. This angle is 
determined through sin^ ^ = Ali^/M"^. Here we do not add the configuration at (/> + vr. The cross 
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section for (/> > vr is obtained by reflection w.r.t. the axis </> = vr. To make sure that we do not get a 
contribution from a region where our perturbative analysis is not reliable we have imposed an abso- 
lute lower cutoff on of 1 GeV^. In addition we have chosen bins in to keep the center of mass 
energy of the photon-pomeron subprocess roughly constant. In Fig. ^ we present the cross section 
with cuts in 9 instead of cuts in imposed. Now the interference term in (6) can no longer be 
neglected, and we expect a slight asymmetry between (p = and (p = ir. In the left diagram we have 
chosen 20GeV^ < M"^ < 50GeV'^ and two different cuts in 9, 7r/16 < 9 < 7r/4 (upper curve) and 
7r/4 < 9 < 7r/2 (lower curve), respectively. Again, as expected, for larger 9, corresponding to larger 
k'^, the asymmetry is more pronounced, but the normalization is smaller. The cos </)-asymmetry 
is clearly visible and it can be seen that the coefficent of the cos (p-teim is positive. In the right 
diagramm we have chosen the same 0-cuts but the M^-bin 50GeV'^ < < 100 GeV^. Here the 
coefficent of the cos(/>-term is negative. This is due to the fact that for large masses the function 
II (eq. dU)) becomes negative. This sign change is a unique property of the interference term. 

Finally we want to illustrate the essential difference between the diffractive cross section based on 
the two gluon exchange model and the boson gluon fusion cross section. We have calculated the total 
cross section for quark-antiquark jet production in ep-scattering based on the photon gluon fusion 
mechanism. The explicit formulae for this cross section can be obtained with the method sketched 
in 3.) and can be found in In Fig. ^ we present the normalized (to unit integral) differential cross 
section in (f) with all other variables integrated using the cuts Q2 > 10 GeV^ ry < 10"^, k^ > 2 GeV^ 
and 50GeV2 < W"^ < 220 GeV^. This curve nicely demonstrates that in the boson gluon fusion 
case the jets prefer a direction in the electron plane (defined by </> = 0), whereas in the two gluon 
exchange case (Figs. ^^) the jets prefer a perpendicular direction. 

5.) In this paper we have calculated the azimuthal dependence of two quark-antiquark jets with 
large transverse momenta in DIS diffractive dissociation. Our main result, contained in eq. (|6|), 
is the striking sign structure of the angle dependent term. It leads to a characteristic maximum 
of the cross section at (/> = 7r/2, quite in contrast to the angular dependence of the photon-gluon 
fusion subprocess which peaks at and vr. An experimental observation of this maximum at 7r/2 
would clearly confirm our present understanding of the 'hard Pomeron'. 



Some time ago it has been suggested that the photon-gluon fusion subprocess, accompanied by 
the additional exchange of soft gluons [ 12 1 or by a suitable modification of the final state interac- 
tions |13| might represent the basic mechanism of the diffractive dissociation observed at HERA. 
In principle, the different dependencies upon (p discussed in this letter might help to discriminate 
between the two mechanisms: the exchange of a color singlet two- gluon state (hard Pomeron) and 
the single gluon exchange in the photon-gluon fusion model. There are, however, two caveats to be 
kept in mind. First, the perturbative calculation of the two-gluon exchange can be justified only 
for jets with large transverse momenta; this excludes any prediction for the diffractive structure 
function which may very well be dominated by soft final states (e.g. Pomeron-remnant jets). The 
discussion in [l^, on the other hand, refers to the diffractive structure function and not to specific 
final states. Secondly, for the case of the photon-gluon fusion we do not know to what extent final 
state hadronization will modify the azimuthal dependence of the partonic cross section: since the 
qq pair is in a color octet state, one expects a color connection with the proton system which may 
affect the angular distribution. In contrast, for the two-gluon exchange model we believe that the 
partonic cross section is rather robust against hadronization effects. Namely, our perturbative cal- 
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culation describes the azimuthal distribution of the qq dipole which initiahy has a small transverse 
size. The subsequent hadronization will be restricted to the qq-system, and thereby not change the 
orientiation. 

In summary, final states in the diffractive dissociation of the photon which contain only hard 
jets represent a novel class of hard processes which are calculable in perturbative QCD: jet pro- 
duction from the 'annihilation' of a photon and a color singlet two-gluon state (a hard Pomeron). 
The qq jets considered in and in the present paper represent the very first step along this line 
- somewhat analogous to the two-jet final states in e^e~ annihilation. What has to come next is 
the order-Os correction to the two-jet configuration, together with the three-jet final state. 
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Figure Captions 

Fig. |l^ : One of the four diagrams contributing to the hard process ^* +p — > qq+P- The outgoing 
(anti)quark momenta are fixed and the same for all four diagrams. The blob depicts the 
unintegrated gluon structure function Tq- 
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Fig. lb : Definition of planes and angles in the 7*iP cms {q + xjpp = 0). The leptonic plane is 
given by the vectors I and q = I — (!.\ the jet plane by the vectors xpp and k). is the angle 
between these two planes. 9 is the polar angle of the jet direction k. 

Fig. |2| : One of the two diagrams for the boson gluon fusion (BGF) process. The square of the 
BGF amplitude again leads to the gluon structure function. 

Fig. ^ : The (j) dependence of the ep - cross section in eq. (^) integrated over the other variables. 
Displayed is the sum of all terms (solid line), and the contribution of the angular independent 
transverse (dotted line) and longitudinal (dashed line) terms. 

Fig. ^ : The <j) dependence of total ep - cross section with different cuts in : 1 GeV^ < < 
2GeV2, 2GeV2 < k^, SGeV^ < k^ (from top to bottom). 

Fig. ^ : The dependence of total ep - cross section with different cuts in 9 and M^. In the left di- 
agram we have 20GeV2 < < SOGeV^ and 7r/16 <9 Ktt/A (upper curve), 7r/4 < 6* < 7r/2 
(lower curve). In the right diagram we have 50GeV^ < < 100 GeV^ and vr/lG < 9 < tt/A 
(upper curve), 7r/4 < 9 < tt/2 (lower curve). 

Fig. ^ : Normalized differential cross section for two jet production based on the boson- gluon 
fusion mechanism. 

Figures 




Figure la: 
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